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ABSTRACT 


Spectra of the archetypal Type II Seyfert galaxy NGC 1068 in a narrow wave¬ 
length interval near 3.7 /zm have revealed a weak absorption feature due to two 
lines of the molecular ion The observed wavelength of the feature corre¬ 
sponds to velocity of -70 km s“^ relative to the systemic velocity of the galaxy, 
implying an outward flow from the nucleus along the line of sight. The absorption 
by along with the previously known broad hydrocarbon absorption at 3.4 yum 
probably are formed in diffuse gas that is in close proximity to the continuum 
source, i.e. within a few tens of parsecs of the central engine. Based on that 
conclusion and the measured Hjj' absorption velocity and with the assumption of 
a spherically symmetric wind we estimate a rate of mass outflow from the AGN 
of ~1 Mq yr“^. 


Subject headings: galaxies: active - galaxies: Seyfert - galaxies: ISM - ISM: molecules 
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Introduction 


Active galactic nuclei (AGNs) contain a number of components on various spatial 
scales: the accretion disk, broad-line region (BLR), torus, narrow-line region (NLR), etc. 
Rather than being the static, well-dehned structures that these terms may imply, in reality 
they are probably complex, dynamic, and in close interaction with each other and with the 
AGN’s host galaxy. For example, spectroscopy of the NLR in some high-redshift quasars 
suggests that an AGN -driven outflow can suppress star formation in the host galaxy 


flGano-Diaz et ah 


2 OI 2 I) . On much smaller scales, the dusty torus of the AGN unihed model 
may be pr oduced by gas flowi ng towards the AGN and play a central role in regulating 


accretion flHopkins et ah 


20121). Gonversely, others have suggested that the torus originates 


in gas flowing out from the AGN (Elitzur & Schlosman 2006), or is formed by the interplay 
between inflow and outflow (Wada 2012). Understanding the motions of gas in AGNs may 
therefore reveal clues about the way in which AGN and their surroundings interact with 
each other and evolve over time. 

One might expect that where there is a dusty torus there are molecules and that 
therefore, for a heavily obscured Type II AGN, infrared spectroscopy of the central source 
would reveal absorption lines from molecules in or associated with the obscuring toroid, 
thereby probing the kinematics of the material on the line of sight to and close to the central 
engine. The most extensive such spectroscopy to date has been of the archetypal Seyfert 


II galaxy, h 


flLutz et ah 


GG 1068. However, tho s e observations, at a variety of spectral resolutions 


2004 


Mason et al. 


2006 


Geballe et al. 


20091) . failed to detect lines of what 


should be the strongest infrared molec ular band, the funda. mental vibration-rotation band 


of carbon monoxide (GO) near 5 /rm. 


Geballe et al. 


(120091) concluded that the gaseous 


environment directly in front of NGG 1068’s AGN, as viewed at mid-infrared wavelengths, 
is more similar to diffuse clouds, which have densities of 10-300 cm“^, than dense clouds. 
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which have higher densities. This conclusion was based not only on the absence of CO, 
but also on the presence of the 3 .4-/ini hydrocarbon absorption feature in NGC 1068 


flBridger. Wright, fc Geballe 


19941) . which is not found in Galactic dense clouds and is a 


reliable signature of a diffuse cloud environment. In Galactic diffuse clouds typically half of 
the hydrogen is in molecular form (H 2 ), but only about one percent of the carbon is in GO. 


A similar percentage of carb on in GO on t 


with its non-detection there flGeballe et ah 


le sig htline to NGG 1068 would be consistent 


2009|). 


Unlike emission spectroscopy, absorption spectroscopy samples material only directly 
on the line of sight to a background continuum source. At 3-5 fim the bright central source 
of NGG 1068 has an unresolved or marginally resolved core of linear dimension less than 
0'.'05 (3.5 pc at the 14.4 Mpc dist ance to NGG 1068) and a complex surrounding morphology 


of characteristic dimension ~1" flGratadonr et ah 


20061) . When observed at 3-5 fim in good 


seeing on large telescopes without the benefit of adaptive optics the nuclear source has 


an effective fu 


(iGeballe et ah 


wid th at half maximum (FWHM) of ~ 0"3 (20 pc) in the N-S direction 


20091 ). The galactic disk of NGG 1068 is inclined at 40° sue. 


side of the minor axis of the galaxy lies to southeast (jde Vaucouleurs et al 


1 that the near 
I 991 I. see also 


de Vaucouleurs 1958); thus the sightline to the AGN passes through only a small portion of 
the galaxy within the nuclear region, and does not pass through the disk of the galaxy. 

Gombined knowledge of the detailed morphology and kinematics of the gas in the 


vicinity of an AGN is importa nt for understanding t 


le in 


central engine. For NGG 1068, 


Garcia-Burillo et al. 


all and/or outflow close to the 


20141) have used the Atacama Large 


Millimeter Array (ALMA) to observe the molecular gas and dust in this region, finding that 
the bulk of each is located in a circumnuclear disk (GND) of outer radius ~200 pc whose 
center is offset by a few tens of parsecs from the position of the AGN. From the profiles 
and intensity ratios of some of the submillimeter lines in the GND they have identihed 
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the kinematic signature of a massive outflow, with a rate of many tens of solar masses per 
year. They identify the gas responsible for this signature as disk material being entrained 
by the ionized wind from the AGN . Closer to the AGN, on scales of 3-30 pc from the 
center, 


Muller Sanchez et ah 


(j2009[) have found that the kinematics and morphology of 


the 2.12 yum H 2 line emission indicate linear streamers of gas are fueling the nucleus, 
and thus presumably the outflow. They interpret some of this gas as being situated on 
th e near side of t h e nuc leus. A contrasting interpretation of the H 2 line emission is given 


by 


Barbosa et af 


(120141) . based on their velocity channel maps in the inner few hundred 


parsecs, obtained at a resolution of 8 pc. They hnd that their H 2 data are consistent with a 
radially expanding ring located in the plane of the galaxy. Regardless of the interpretation 
it is likely that the H 2 line emission arises in a small fraction of the molecular gas that is 
at high temperature immediately following shock-excitation or perhaps X-ray heating, and 
whose motions are not necessarily representative of the kinematics of the neutral nuclear 
gas as a whole. 


The CND observed by 


Garcfa-Burillo et ah 


(120141 ) overlaps the AGN on the sky, so it is 


important to know whether it passes in front of, behind, or through the AGN. The nucleus 
of NGC 1068 is distinguished by a prominent bicone of ionized gas at a position angle 
of ~10°, extending several arc seconds on either side of the center. The ionization cone 
extending NNE from the AGN is much more prominent at optical wavelengths than the 
southern cone. Both that and the blue-shifts of the emission lines observed in it indicates 
that the NNE cone is tilted toward the sun. The plane of the GND is expected to be roughly 
perpendicular to the bicone. If so then the GND lies largely or entirely behind the AGN 
where they overlap on the sky. Nevertheless, in view of the huge quantity of molecular gas 
and dust that are present in the GND as well as the asymmetric location of the GND with 
respect to the central engine it remains surprising that until now no molecular absorption 
lines had been seen on the line of sight to the AGN. 
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Clearly, an observation of absorbing gas on the line of sight to the AGN could provide 
important additional information about gas motions close to the central engine. Because 
the search for absorption lines of CO has been unsuccessful, an alternative infrared probe 
of the foreground gas column is needed. One possible candidate species is Hg . Although it 
is obviously an ion, Hjj’ is created by the reaction H 2 + Hj —> + H and resides in 


largely neutral and at least partially molecular gas. Vibration-ro tation lines of I 


generally weak, are readily detectable no 


only in dense clouds fiGeballe fc Oka 


, alt hough 


19961) . but 


also in numerous Galactic diffuse clouds 

McGall et ah 

2002 ; 

Indriolo et ah 

2007 

), including 

e.g. Gygnus OB2 No. 12 ( 

Geballe et ah 

1999 

), whose 3.-^ 

:-wm feature is less than half the 


iwhittet et ah 

1997; 

Geballe et ah 

2009) 


20091) . Unless the 


lines of Hg are greatly broadened toward NGC 1068 one would expect to detect them on 
the line of sight to the AGN. The close pair of vibration-rotation lines (1-0 i?(l,l)“ and 1-0 
i?(l,0)), at rest wavelengths of 3.66808 pm and 3.66852 pm, respectively, which arise from 
the lowest levels of para and otho Hg", are the best candidates for detection. In both diffuse 
and dense clouds in the Galactic plane these are the only two levels that are appreciably 
populated and the doublet is the strongest Ht absorption feature. 


2. Observations and Data Reduction 

L-band spectra of NGC 1068 were obtained at the Frederick C. Gillett Gemini 
North Telescope on UT 2011 July 24 and October 11 and on UT 2014 November 28 and 
December 3, for programs GN-2011B-Q-53 and GN-2014B-Q-57, respectively. Both sets 
of observations used the Gemini North Infrared Spectrograph (GNIRS), but in different 
conhgurations. The 2011 spectra were obtained using the 111 1/mm grating, short focal 
length camera and the 0.45 " wide slit and covered 3.49-3.79 pm. The 2014 spectra were 
obtained through the same slit, but used the long focal length camera and the 111 1/mm 
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grating, which covered 3.64-3.74 and were heavily oversampled in wavelength. In 
each case the resolution at the wavelength of the sought after lines was 0.00082 yum 
(70 km s“^). The observations were obtained in the standard stare/nod-along-slit mode, 
with a nod of -|-/- 3". The slit was oriented along the axis of the bicone at position angle 
15° (E of N). The total exposure times were 64 minutes in 2011 and 120 minutes in 2014. 
Early type telluric standard stars were observed at closely similar air masses to NGC 1068 
either before or after each set of observations. The measurements were made in photometric 
or near photometric conditions and in good seeing. 

The initial stage of data reduction included flat-fielding, spike removal, extraction 
of spectra over 0'.'75 long regions of the slit, wavelength calibration (using telluric lines 
and accurate to 0.0001 yum, corresponding to 10 km s“^ ), and ratioing by the spectra of 
the telluric standards. Normalized ratioed spectra were then shifted to compensate for 
differences in the earth’s orbital motion relative to the coordinates of NGC 1068 on the 
different observing dates, and then combined to produce one spectrum for each of the two 
programs. The continuum of each spectrum had slight curvature, which we suspect is an 
instrumental effect due to differences in the distribution of light from NGC 1068 and the 
standard stars and different positioning of them in the slit; to remove it and obtain flat 
spectra each spectrum was divided by a spline fit to its continuum. The resultant spectra 
were then resampled on a common wavelength scale and coadded with equal weighting (as 
their noise levels are similar) in the narrow wavelength interval encompassing the expected 
wavelength of the Hj]' doublet, near 3.682 yum. Strong telluric absorption lines are present 
on either side of this wavelength, but the interval 3.675-3.686 ^m, where the redshifted Hjj' 
doublet is expected to fall, is devoid of significant interference. The spectra from semesters 
2011B and 2014B and the final combined spectrum are shown in Eig. 1. 
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Wavelength (vae. microns) 


Fig. 1.— Spectrum of the infrared nuclear source in NGC 1068 near the Hjj" para and ortho 
line pair, at a resolution of 70 km s“^. The spectra from 2011 and 2014 (light dotted lines) 
and the combined spectrum (dark line) are shown, as is a 100 times scaled down spectrum 
of the atmospheric transmission (dark dotted line).The wavelengths of the Hg" lines and of 
the adjacent [A1 vi] coronal line are indicated for a systemic heliocentric velocity of 1142 km 
s“^ (see text). 
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Results 


The combined spectrum in Fig. 1 shows a detection of a very weak absorption feature 
with its centroid at 3.6815 /rm, which we identify as the blended pair of lines of para 
an d ortho Hj~ m e ntione d above. A coronal emission line of [Alvi], previously reported 


by 


Geballe et ah 


(120091) is also present in the spectrum. Based on the fluctuations in the 


adjacent continuum in the combined spectrum the detection of the Hj]' feature is at the 
6(7 conhdence level. However, the differences in the 2011 and 2014 spectra on the long 
wavelength portion of the continuum suggest that the true conhdence level is not quite 
that high. In estimating the continuum and the noise (the latter based on point-to-point 
huctuations in the spectrum) we ignored the spectral intervals near the telluric lines at 
3.674 pm and 3.688 pm, as the noise levels there (which are dominated by huctuations in 
the sky background) are signihcantly higher than elsewhere. 

Assuming a systemic heliocentric velocity of 1142 km s“^ (from the H I Parkes All 
Sky Survey, as reported in the NASA/IPAC Extragalactic Database), the centroid of 
the absorption is blue-shifted by approximately 70 km s“^ from the systemic velocity of 
NGC 1068, indicating that the gas associated with the is howing outward from the 
nucleus. The feature has a measured equivalent width of 5.2 ± 0.8 x 10“® pm and an 
observed FWHM of ~150 km s“^. Assuming that the two lines, which are separated by 
36 km s“^, contribute equally to the feature, as they approximately do in Galactic diffuse 
clouds, the combined column density in the (1,1) and (1,0) levels is 7 x 10^^ cm“^ and the 
intrinsic FWHM of a single absorption line is approximately 100 km s“^. T hat width is 
consid erably narrower than the more highly blue-shifted [Alvi] line in Fig. 1. 


Geballe et ah 


(120091 ) found that the peak intensity of the infrared coronal lines occurs a few tenths of 
an arc-second to the north of the infrared continuum peak. The present data also show 


that separation. Glearly the gas responsible for the absorption by H^ and the gas i 


m 
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which the line emission by [A1 vi] arises are in completely different locations and physical 
environments. 

Other key lines of that could confirm the identification of the 3.6815-yum feature as 
due to and could provide constraints on the temperature and density of the material 
in which it originates are the and i?(3,3)^ lines at 3.71548 /xm and 3.53366 pm, 

respectively. Unfortunately each of these lines is likely to be weaker than the unresolved 
doublet and each is shifted into a wavelength interval where strong telluric absorption lines 
are present. Despite the lack of conhrmatory evidence we are conhdent in the identihcation 
of H 3 +, as ( 1 ) there are no other likely identihcations for the feature and ( 2 ) the line pair 
was expected to be readily detected. For example, if the equivalent width of the doublet in 
NGC 1068 scaled as the ratio of equivalent widths of the observed 3.4 pm and doublet in 
the Galactic center (see Fig. 2), the equivalent width of the doublet in NGG 1068 would be 
3 X 10“^ pm, six times greater than that observed. The surprising weakness of the doublet 
in NGG 1068 is discussed below. 


Discussion 


4.1. Extinctions to and Locations of Continuum and Absorption Feature 


L-band spectra of NGG 1068 have yielded two absorption features, a prominent 3.4 pm 
feature generally ascribed to aliphatic hydrocarbons, and a weak absorption due to a pair 
of closely spaced lines of Hjj". Understanding where on the line of sight these absorptions 


occur requires first of all knowledge of the 
continuum. The adaptive optics images of 


ocation of the source or sources of the infrared 


Gratadour et ah 


( 120061 ) show that within the 


portio ns of the slit from which the spectra of these features were extracted, by 


Geballe et ah 


f 2 OO 9 I) for 3.4 pm feature (0'.'2 x I'.'O) and by us for the H 3 feature (0'.'45 x 0'.'75), the 
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continuum is produced by a high surface brightness core and a complex morphology of 
more extended and lower surface brightness regions with the two being roughly comparable 
in total signal. 


Whatever the explanation for the L-band morphology; the likeli hood is that the 


conti nuum is thermal emission from clouds containing hot dust (e.g., 


Gratadour et ah 


2006h . presumably heated by the central engine. If so, the continuum probably arises 


m 


the vicinity of each r=l (in the L band) surface in this complex region; in other words at 
distances of a few parsecs to a few tens of parsecs from the engine. One can also estimate 
the mean optical depth to the L-band continuum from the depth of the 3.4 /rm feature, 
assuming the relationsh ip between the dept h of the feature and the optical extinction in 


Galactic diffuse clouds flGodard et al 


of the feature is 0.08 flGeballe et al 


2012h holds for NGG 1068, where the optical depth 


20091 ). This assumption yields a visual extinction 
of ~21 mag, which, for normal extinction curves (e.g., extinction proportional to A“^'^) 
implies an L-band optical depth close to unity. The near equality of the optical depths 
derived in these different ways is either a coincidence or, as we suspect, an indication that 
the dust producing the L-band continuum and the carrier of the 3.4 /rm absorption (and by 
association the absorbing H^) are physically close to one another, placing the absorbing 
material within a few tens of parsecs of the central engine. This proximity is not surprising 
in view of the orientation of the disk of the galaxy, as discussed earlier. 


4.2. Low Hg Equivalent Width 

Figure 2 is a plot of the equivalent width of that portion of the Hj]' line pair arising in 
diffuse gas versus the equivalent width of the 3.4-/im absorption feature, for four sources: 
the Galactic plane source Gyg OB2 No. 12, the Galactic center’s tight cluster of four 
sources known as GGS 3 (of which the brightest component is GGS 3-2), the ultraluminous 
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1 10 
W^(3.4jum) (10~^juTn) 


Fig. 2.— Equivalent width of that portion of the absorption by the para-ortho line pair 
ascribed to diffuse gas vs. equivalent width of the 3.4 pm feature, for Cygnus OB2 No. 12, 
the Galactic center Quintuplet source(s) GCS 3, the ultraluminous infrared galaxy (ULIRG) 
IRAS 08572-1-3915, and NGG 1068. Error bars are ±lcr. The diagonal line is a one-to-one 
relationship between the two equivalent widths. The equivalent width of the feature in 
GGS 3 is based on the measurement toward GGS 3-2. In all but GGS3-2 the absorption by 
is assumed to be formed entirely in diffuse gas. In GGS 3-2 the contribution by Hjj" in 
dense clouds has been subtracted. See text for details. 
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infrared galaxy (ULIRG) IRAS 08572+3915, and NGC 10 68. The Hj~ eq u ivalen t widths for 


the fir st thr ee sources are obtain ed from data or hgures in 


Geballe et al 


(120051), and 


Geballe et al 


were m easu red from spectra in 


fll999h . 


Oka et al. 


(12006), respectively. 


Sandford et ah 


he equivalent widt hsoft h e 3.4-Mm featu re 


(120061), and 


Mason et ah 


19911) ■ 


Ghiar et ah 


( 120001 ) ■ 


Geballe et ah 


(120061 ) ■ respectively. For Gyg OB2 No. 12, IRAS 08572+3915, and 


NGG 1068 either all of 


gas; (see 


Geballe et al. 


he H j~ or the vast majority of it is known or thought to be in diffuse 


2006 


for the argument pertaining to IRAS 08572+3915). In the 


case of GGS 3-2, how ever, the sightline passes through both diffuse and dense interstellar 


gas. 


Oka et al. 


(120051 ) have determined the column densities in each environment. A further 


complicat ion is that t 


However, 


Oka et al. 


le lin es are Doppler broadened so that the pair is observed as a blend. 


allowing the equivalent width of the b 
to be determined using the formula in 
to line equivalent width. 


(1200511 also estimated the populations in the (1,1) and (1,0) levels, 
ended pair of lines, w hich absorb from those levels. 


Geballe fc Okal (119961) relating level column density 


If the bulk of the Hg is in the lowest para and lowest ortho levels, as is generally the 
case, then the equivalent width of the line pair represents the total column density of H^. 
One then would naively expect that the equivalent widths in Fig. 2 would scale with one 
another. NGG 1068 clearly deviates the most from following that trend; the discrepancy 
is roughly a factor of six. This suggests either that an unusually low fraction of its total 
Hd" column density resides in its two lowest energy levels or that it has a highly enhanced 
abundance of the carrier of the 3.4-/rm absorption feature, or both. There are several 
possible explanations for this anomaly, which we now list and briefly discuss. 


1. in other energy levels: For cold diffuse gas such as that in front of Gygnus OB2 No. 
12 the (1,1) and (1,0) levels are the only ones that are signihcantly populated. In the 
Galactic center, where the density of the diffuse gas containing the Hjj’ is ~100 cm“^ 
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and the gas temperature is ~250 K the metastable (3,3) level contains about 23% of 
altogether all o f the levels higher than (1,1) and (1,0) contain about 30% 


the Hjj’ and 
of the Ho 


Oka et ah 


2005|) . Thus the correction, which would move its data point 


upward by 30% in Fig. 2 is modest. In IRAS 08572+3915 the large column densities 
implied by the strong line s suggest that the H j" is distributed along a column 


length of hundreds of parsecs (IGeballe et ah 


20061) . Much of this extended region is 


likely to be at temperatures that are no higher than those in the Galactic center, and 
thus we believe that the upward correction is also small. 

:ed with the 


If the Hg observed toward NGG 1068 is in gas physically associa 
contin uum source, as we suspect, its temperature could be higher. 


Gratadour et ah 


(120061) suggest dust te mperatures of ~600 K for the dust emitting in the L band. 


Alloin &: Marcol (119971) . observing at 77, L, and M, found a mean dust temperature 


of 750 K for the central ~l-arc second region. 


Jaffe et ah 


10 nm, hnd a 320 K in the central few parsecs, as do others (iMason et ah 


(I 2 OO 4 J). observing near 


Poncelet. Sol, fc Perrin 


2006 : 


20081) . Thus a signihcant portion of the absorbing column may 


have a higher proportion of in higher levels than the Galactic center. However, 
due to the wide spacing of rotation energy levels in the ground vibrational state 
(e.g., the J=4 levels have excitation energies >1,000 K) as well as depopulation via 
spontaneou s decay of many o f the excited rotational states at the low densities of 


diffuse gas ( Oka Sz Ep 


20041) . we expect the correction to be less than a factor of 


two. Thus we do not think that large populations in higher energy levels can explain 
the discrepant data point in Fig. 2. 

2. Higher gas density and shorter column length: In diffuse gas the density of Hjj' does 
not track the total number density of the gas in which it is located, but is a constant 


whose value depends on the rate of cosmic ra y ionization o 


dissociative recombination of H^ on electrons (IGeballe et ah 


H? a nd the rate of 


1999 ). Thus the column 
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density of scales with column length rather than with gas column density. If the 
mean density of the diffuse environment in front of the AGN in NGC 1068 is higher 
than in typical diffuse clouds, the column density of there will be lower relative 
to other species. Likewise, a larger than normal density of free electrons, which is the 
principal destroyer of Hjj' in diffuse gas, would result in a lower steady state density 
of and also in a lower column density. 


3. Low fractional abundance of H 2 : As discussed in Section 2.3.2 of 


Oka 


(1201311 the 


abundance of Hjj" is approximately proportional to the square of the fractional 
abundance of molecular hydrogen, /(H 2 ) = 2?7,(H2)/[n(H) + 2n(H2)]. A value of 
/(H 2 ) that is smaller by a factor of 2.5 than in typical Galactic diffuse clouds could 
explain the low abundance of H^. In NGG 1068 the abundance of H 2 might well be 
lower than normal in the hostile environment close to the AGN. 


4. High abundance of carbon and/or the carrier of the 3.4-yum feature: 

The 3.4-//m feature is generally attributed to G-H stretching vibrations in aliphatic 
hydrocarbons. As such the strength of the feature is sensitive to the carbon abundance 
and is likely to be greater toward the centers of galaxies, paralleling expected or 
observed increases in metallicity. In diffuse clouds, where most of the carbon is singly 
ionized, a higher abundance of carbon also results in a higher electron density and a 
lower steady state abundance of Hg . Judging from Fig. 2, if these effects alone were 
to account for the discrepancy the carbon abundance in the nucleus of NGG 1068 
would need to be several times higher than in the nuclear regions of IRAS 08572+3915 
or the Galactic center. We have found no observations of NGG 1068 in the literature 
that are pertinent to this issue. 


To summarize this subsection, there are several possible explanations for the unusually 
low value of the ratio of equivalent widths of the and 3.4-/im absorption features on the 





sightline to the AGN in NGC 1068. One, some, or all of them could contribute and we are 
unable to select a likely dominant contributor. 


4.3. Outflow and Mass Loss 

The observation of a blue-shifted Hjj" absorption feature at 3.7 ^m indicates that on 
the line of sight to the nucleus, and thus quite far removed from the axis of the ionization 
bicone, there exists outflowing gas located close to where L-band continuum is formed, at a 
distance of possibly only a few parsecs or a few tens of parsecs from the central engine. The 
relationship of this outflow to the apparently much more powerful flow largely within the 
ionization bicone, which probably originates much closer to the central engine, is unclear. 

Based on the absorption feature, a very crude estimate of the rate of mass outflow 
can be derived. We based our estimate on (1) the Hg being in the same diffuse interstellar 
environment that produces the 3.4-/im hydrocarbon feature and (2) our interpretation that 
the similar values of the likely optical depth to the source of the L-band continuum and 
the optical depth estimated from the hydrocarbon feature mean that the diffuse gas is in 
proximity to the continuum sources rather than at a great distance from them. Using a 
diffuse cloud gas density of 100 cm“^ and a mean distance of the flow at that density of 
20 pc, together with the highly dubious assumption of isotropy of the outflow, we obtain 
an outflow rate of 1 Mq yr“^. The rate is proportional to the square of the assumed 
mean distance, which together with the assumption of spherical symmetry have high 
uncertainties. Other contributors to the uncertainty are the outflow velocity, which depends 
on the choice of systemic velocity, and the mean density of the absorbing diffuse gas. We 
regard the above value of mass outlaw rate as an order of magnitude estimate only. 

If our estimate of 20 pc for the location of the diffuse gas is roughly correct then the 
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outflow observed in Hj]" exists on distance scales smaller than those probed by the ALMA 
millimeter and submillimeter spectroscopy of the CND . Based on the analy s is of m olecular 


line prohles at distances of 100-200 pc from the center 


Garcfa-Burillo et al.l f)2014j) recently 


deduced the existence of an outflow originating in the AGN with a mass loss rate 1-2 
orders of magnitude greater than we derive from the blueshifted Hjj' absorption. They also 
mea sure outflow ve l ocities of up to a few hundred km s“^, considerably higher than seen in 




Barbosa et ah 


(120141) also observe higher velocities, but derive a mass loss rate, based 


on the kinematics of the ionized gas similar to our rough estimate for the diffuse molecular 
gas. Unlike the outward flowing observed by us, which is far from being along the axis of 
the bicone, the outflowing gas that both of these groups have observed is within or at least 
much closer to the bicone. Much higher velocities, up to 1,000 km s ~^, also are observed 


from the bicone in numerous optical and infrared em ission lines (e.g.. 


Grenshaw fc Kraemer 


2000 ; 

Das et ah 

2006; 

Poncelet. Sol. & Perrin 

2008) 


The differences in mass loss rates derived here and by 


Garcfa-Burillo et ah ( 2014 ) 


thus may indicate large anisotropies in the outflow, but also could easily be due to time 
variability both in intensity and in the direction of the outflow. Approximately 10^ years 
are required for gas moving at 1,000 km s“^ to travel 100 pc (the distance from the 
central engine to the inner edge of the GND). On the other hand, it seems reasonable that 
although the outflow is primarily within the bicone, it has an associated and much weaker 
component in other directions. That portion of the flow would be more shielded from the 
strong ultraviolet held present in the bicone and would allow gas that on the line of sight 
to become or to remain partially molecular, although quite close to the AGN. 
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Conclusions 


We have detected a weak absorption by blue-shifted Hg on the line of sight to the 
central bright infrared continuum source of the type II Seyfert galaxy NGC 1068. This is 
only the second detection of in an extragalactic object. Surprisingly is the only 
interstellar molecule in NGC 1068 that has been detected in absorption toward its AGN. Its 
presence is consistent with the existence in front of the AGN of a significant column density 
of diffuse gas, previously known to be present by a prominent 3.4-yum hydrocarbon feature, 
and also from the absence of absorption in the fundamental band of carbon monoxide. 

We tentatively conclude that the diffuse gas is situated within a few tens of parsecs of 
the central engine and estimate an (isotropic) mass outflow rate of 1 Mq yr“^ with an 
uncertainty of an order of m agnitude. This rate i s lower by 1-2 orders of magnitude than 


the rate recently deduced by 


Garcfa-Burillo et al.l (120141) from ALMA observations of line 


pr ohles in the 100-200 p c radius CND, but is roughly consistent with the rate derived 


by 


Barbosa et al 


(12014^ . The Hg absorption is considerably weaker than expected from 


the strength of the 3.4-yum feature. This could be explained by the gas being at a higher 
temperature, by the fraction of hydrogen in molecular form being 2-3 times lower than 
typical for diffuse clouds, by an unusually high abundance of carbon in the nucleus, or by 
the diffuse gas in which the absorption takes place being considerably higher density than 
standard Galactic diffuse clouds in which Hg has previously been detected. 
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